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Abstract 



We present the calculation of the two-loop spin splitting functions {x) = q,g) con- 
tributing to the next-to-leading order corrected spin structure function gi{x,Q'^). These 
splitting functions, which are presented in the MS scheme, are derived from the order 
contribution to the anomalous dimensions 7™ = q,g). The latter correspond to the local 
operators which appear in the operator product expansion of two electromagnetic currents. 
Some of the properties of the anomalous dimensions will be discussed. In particular we find 
that in order the supersymmetric relation jVl + 7™ — 7™ — 7™ = is violated. 



1 Introduction 



During the last few years there has been a great deal of activity in the area of polarized 
lepton-hadron physics both from the experimental as well as the theoretical side. This interest 
started with the discovery of the EMC-experiment Q that the Ellis- Jaffe smn rule @, which 
represents the first moment of the spin structure function gi{x,Q'^), was violated by the 
combined SLAC-EMC data |^. This discrepancy between theory and experiment, also 
called the "spin crisis" , came as a great surprise because one expected that sum rules derived 
in the context of the constituent quark model, which is valid at low energy scales, should also 
hold at large energy scales characteristic of the current quark (parton) regime. In particular 
the constituent quark model assumes that the spin of the proton can be mainly attributed 
to its valence quarks and the sea quark contribution is negligible small. This assumption 
leads to a value of the Ellis-Jaffe sum rule which is appreciably larger than the one found by 
experiment. Although more recent experiments @, ^ lead to a result which is closer to the 
theoretical prediction the discrepancy is still large enough to warrant explanation. 

Many theorists have tried to explain the above discrepancy (for recent reviews see Q) in 
the framework of perturbative and also non-perturbative QCD. From this theoretical work 
one can draw the conclusion that the interpretation of the spin structure function gi{x,Q'^), 
using the ideas of the operator product expansion (OPE) and the QCD improved parton 
model, is not as simple as that given to the structure functions which show up in unpolar- 
ized lepton-hadron scattering. In particular the axial vector operator is renormalized due 
to the Adler-Bell-Jackiw anomaly. Therefore the interpretation that the polarized parton 
densities represent the spin carried by the corresponding partons does not hold anymore. 
Fortunately this operator cancels in the Bjorken sum rule Q so that the latter has a more 
reliable theoretical basis. It is therefore no surprise that its result is in agreement with recent 
data i, |,|]. The ab ove theoretical work also led to many different parametrizations of the 
parton densities in terms of which the spin structure function gi{x,Q'^) can be expressed. 
One of the key issues is the role of the gluon density which can account for the negative 
contribution to the Ellis-Jaffe sum rule depending on the chosen scheme. However if one 
wants to give a complete next-to-leading order (NLO) description of gi{x,Q'^), and not only 
its first moment, one needs a full knowledge of the order coefficient functions, which are 
known (see e. g. |lO|, |ll|) and the order corrected Altarelli-Parisi (AP) spin splitting 
functions Pij = q,g). The lowest order AP-splitting functions pj^j^ have been calculated 



in ||12[ and [|13[ respectively using different methods. In the operator product expansion 
(OPE) techniques are applied to obtain the anomalous dimensions of the composite opera- 
tors appearing in the spin dependent part of the current-current correlation function. The 
latter appears in the expression for the deep inelastic cross section. The authors in [13| 
have used the parton model approach. The NLO (order a^) splitting functions Pql^'^ and 
Pql^ have been computed in Q using the standard techniques of perturbative QCD. They 
emerge while performing mass factorization on the order corrected parton cross sections 
of the processes 7*5 and j*g which contribute to the deep inelastic spin structure function. 
Unfortunately the remaining splitting functions Pgl^ and Pgg^ could not be obtained in this 
way since they do not show up in the mass factorization of order corrected parton cross 
sections. This can be traced back to the phenomenon that there is no direct coupling of 
the virtual photon 7* or any other electroweak vector boson to the gluon. Therefore Pgl^ 
and Pgl^ will appear in the mass factorizaton of the order corrected parton cross sections 
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which are very difficult to calculate. In order to avoid the above complication we will resort 
to the standard OPE techniques to calculate the missing splitting functions which are de- 
rived from the inverse Mellin transform of the anomalous dimensions of composite operators. 



The paper is organized as follows. In section 2 we introduce our notations and present a 
short discussion of the composite twist-2 operators contributing to the spin structure function 
gi{x,Q'^). Here we also derive the general form of the renormalized and unrenormalized 
operator matrix elements (OME) where the operators are sandwiched between polarized 
quark and gluon states. The calculation of the OME's is presented in section 3, from which 
one extracts the anomalous dimensions and the AP splitting functions which are presented 
in the MS scheme. Further we give the lowest order coefficient functions of gi{x,Q'^) in the 
same scheme. The properties of the anomalous dimensions are discussed in section 4. In 
Appendix A one can find the operator vertices needed for the computation of the operator 
matrix elements in section 3. The tensorial reduction of the Feynman integrals which show 
up in the calculation is discussed in Appendix B. 

2 Operators contributing to the spin structure function 

gi{x,Q'^) 

In this section we specify the composite operators which appear in the light-cone expansion of 
two electromagnetic currents. Furthermore we present the operator matrix elements (OME's) 
as a power series in the strong coupling constant. The coefficients of the perturbation series 
are determined by the renormalization group (Callan-Symanzik) equations. We will write the 
OME's in the most general way so that they can be used to extract the anomalous dimensions 
of the composite operators. The light-cone expansion of two electromagnetic currents is given 
in 1 12 1 and reads as follows 
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z,^Rr-'-{Q). (2.1) 



In the above we only consider the contribution of twist-2 operators. The index i of the locally 
gauge invariant operators Q^^ and ^^('^■■•m™ stands for the representation of the flavour 
group SU{nf). Notice that the operators are also irreducible representations of the Lorentz 
group which means that they are traceless and symmetric in the Lorentz indices /ii • • -/im. 
The Wilson coefficient functions, denoted by C^™ {k = 1,2) and £'j"j, can be expressed into 
a perturbation series in the gauge (strong) coupling constant g. Notice that all the above 
quantities are renormalized which is indicated by the renormalization scale /i. The product 
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of the two electromagnetic currents appear in the hadronic tensor defined in polarized deep 
inelastic lepton-hadron scattering which is given by 

W^,{p,q,s) = ^ f d^ze"i^p,s\J^{z)MO)\p,s) 



47r 

= W^,ip,q)+iW^,ip,q,s). (2.2) 

Here p and s denote the momentum and spin of the hadron respectively and q stands for the 
virtual photon momentum. The symmetric part of the hadronic tensor is given by 

W^AP^Q) = [-9i^u + -^)Fi{x,Q ) + [Pf, jqf,)[Pu jl") ' ^•^) 

^ q"^ q'^ q'^ pq 

while the antisymmetric part is equal to 



W^^{p,q,s) = -—Cf^uapq 



(2.4) 



pq 

with the properties s ■ p = 0, = 1 and m denotes the mass of the hadron. The Bjorken 
scaling variable is given by X = Q^/{2pq) and = -q^ > 0. The spin averaged structure 
functions are denoted by Fi^{x,Q'^){k = 1,2). In polarized electroproduction one has in 
addition the longitudinal spin structure function g\{x,Q'^) and the transverse spin structure 
function g2{x-,Q'^)- The twist-2 operators O^^"''^™'(0) corresponding to the spin averaged 
structure functions are given in the literature and their anomalous dimensions have been 
calculated up to two-loop order p^ j-j]!?!. The twist-2 operators contributing to the spin 
structure functions are given by |jl2| 

= {(^(^)757'^^i^'^^---^'^'"^A,V'(^)-(traces)|, (2.5) 

^Mr-/^m(^) = i« 5 {(V; (2)757/^1 . . _ (traces)} , (2.6) 

R^s]g^^{z) = i"" S ^^^e'''"^^Tr{Ffs^{z)D^'^ ■ ■ ■ D"^-^ F^^ (z)) - (traces)| . (2.7) 

The symbol S in front of the curly brackets stands for the symmetrization of the indices 
Hi ■ ■ ■ Hm and Aj is the flavour group generator of SU{nf). The quark and the gluon field 
tensor are given by tpiz) and F^^,{z) respectively and F^^i, = FI^^^T"" where stands for 
the generator of the colour group SU{N) {N = 3). The covariant derivative is given by 
F>ij, = d/j, + igT"" A'^{x) where A'^^[z) denotes the gluon field. From eqs. ( |2.5D -(^^) one infers 
that with respect to the flavour group one can distinguish the local operators in a non-singlet 
part represented by RNS,q and in a singlet part consisting of Rs,q and Rs,g- 

In the Bjorken limit (Q^ 00, x = Q'^/{2pq) flxed) the current-current correlation func- 
tion in (^]^) is dominated by the light cone = so that it is justified to make a light cone 
expansion for the product of the two electromagnetic currents. When Q'^ — > 00 the lead- 
ing contribution of gi{x,Q'^) consists of the twist-2 operators listed in (^)- (l2.7D whereas 
5(2 (x, Q^) also receives contributions of twist-3 operators which are not given in the expansion 
in eq. ( |2.lD . Since we are only interested in the longitudinal spin structure function gi{x, Q^) 



we can limit ourselves to the renormalization of the twist-2 operators mentioned in (2.5)- 
(^). Inserting the light cone expansion for J^{z)Jy{d) in ( |2.2| ) one can derive the following 
relation 

[\xx"'~^gl{x,Q'')=Y^AT{p\^Ji^g)E^l{Q^^^\g) ; m odd. (2.8) 

JO : 
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The left-hand side of the above equation stands for the Melhn transform of gi{x, Q"^) and the 
right-hand side is given by the operator matrix element (OME). 



{p,s\R 



{0)\p,s) =i'"Ar(P •••P^™) - (traces)} 



with i =NS,S, and Ef^ stands for the coefficient function. 



2\m 



Q \ m— 1 



d ze 



iqz 



lezQ 



E^^iz^ -iezo,fi^g). (2.10) 



The Q'^-evolution of the spin structure function is determined by the anomalous dimensions 
of the composite operators in eqs. ( |2.5[) - (|2.7| ). They are obtained from the renormalized 
partonic OME's 

{j,P,s\R^,y-^-\j,P,s)=A^^,y,f,^g)S{{s'^^p^'---p^-^)-{tv^^^^^ (2.11) 

where now the quark and gluon operators are sandwiched between quark and gluon states. 
The indices in ( 2.11 ) stand for k = NS, S and i = q, g; j = q,g- The A"/^- are derived from 
the Fourier transform into momentum space of the connected Green's functions 



(O|r(<^,(x)<V'^-(O)0,(y)|O), 



(2.12) 



where the external lines are amputated. The fields 0j(x) stand either for the quark fields tp(x) 
or for the gluon fields A'^(x). The renormalized partonic OME's satisfy the Callan-Symanzik 
equations 



^^4- + Pi9)^ + ^ia,g)^+-f'^s,qq(■9) 



9/i 



da 



^Sv5,gg(f'^/^^5,a) = 0, 



and 



(m|; + /3(5)^ + 6{a,g)-^)6i, + j^^.^ig) 



4^.i(p^^^9,a)=0. 



(2.13) 



(2.14) 



Here l3{g) denotes the /3-function which in QCD is given by the following series expansion 



Pig) 



9 



167r2 



Pi 



9 



+ ... 



11 4 
Po = —CA--Tfnf, 



Pi 



34 



20 



^Ci-ACFTfnf--CATfnf. 



(2.15) 

(2.16) 
(2.17) 



Further 5{a^g) is the renormalization group function which controls the variation of the 
OME's under the gauge constant a. Choosing the general covariant gauge one obtains in 
QCD the following result 



5{a,g) 



-az. 
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"16^2 



where Za is given by 
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+ 



(1 - a) ) + -Tj. 



(2.18) 
(2.19) 
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Furthermore the colour factors of SU(N) are defined by Ca = Cp = {N"^ — 1)/{2N), 
Tf = 1/2 and nj stands for the number of hght flavours. The anomalous dimensions are 
given by the series expansion 

2 

I +■•• (2.20) 



Ik. 



(0),m 9^ (l),m ( JJ^ 



167r2 



Notice that for the subsequent part of this paper we do not need higher order terms in 
P{g), 6{a,g) and jJP'^j- As an alternative to using the renormalized OME's the anomalous 
dimensions can also be derived from the operator renormalization constants ZJT'-j which relate 
the bare operators Ri^k to the renormalized operators Ri^k Q- The renormalization of the non- 
singlet operator proceeds as 



z 



(2.211 



Since the singlet operators in ( |2.6D and (2.7) mix among each other the operator renormal- 
ization constant becomes a matrix and we have 



Zf,Ae,g)R''s'r-iz). 



's,ij\^^yj^''S,j 

Now the anomalous dimensions also can be obtained from 



(2.22) 



where 



d 



(3ig,e) {Z^% 



1 



dZsj 



dg 



(5{g,e) = -eg + P{g) 



(2.23) 



(2.24) 



Here e = n — 4 indicates that we will use n-dimensional regularization to regularize the 
ultraviolet singularities occurring in Zk^ij which are represented by pole terms of the type l/e^. 
The computation of the OME's proceeds in the following way. First one adds the operators 
(^)-(|2y7|) to the QCD effective lagrangian by multiplying them by sources Jfii---^im{z)- The 
calculation simplifies considerably if the sources are chosen to be equal to J, 



A^, • • • A^^ with A2 



Ml ■■■Mm " 



0. In this way one eliminates the trace terms on the right-hand side of 
eq. ( |2.11| ). The Feynman rules for the quark and gluon operator vertices are given in Appendix 
A. Starting from the bare lagrangian, which is expressed in the bare coupling constant and 
bare fields and operators, one obtains the following general form for the unrenormalized 
OME's. For the non-singlet OME we have 

e/2 



A 



NS,qq 




„(o) 



TAfS,-?? p + ^N's,qq + ^ ^NS,qq 



n (0) 



+ Tns, 



(0) „(i) 



da 



qq'^NSm 



a 



(1) 

NSm 



da 



(2.25) 



^Notice that in the subsequent part of the paper the unrenormalized quantities will be indicated by a hat. 
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Here 



Ss is a factor which originates from n-dimensional regularization. It is defined by 
The singlet quark OME can be written 



as 



As,qq = ANS,qq + Aps,qq, 

where the pure singlet (PS) 



ApS,qq = 



9 



part is given by 

,2 \2 



167r2 



S 



_p2 
.7^ 



'1 (0) (0) J_ 
2'yS,qg^S,gq ^2 



+ I 2^PS,qq + ^S,qg(^S,gqj ^ + «P5,g? 



SO that 7gjg = 7^^^^^ + jPsm- 
The other OME's can be expressed in the renormalization group coefficients as 



^S,qg = 



+ 



+ 



^167r2 



1^167r2 



1 (0) / (0) , (0) \ a 
^Ts/qg [rS,qq + rS,g9)-^^ 

1 (1) na (1) I (0) (1) I (0) (1) ' 



^(0) 1 1 

IS,qg J 



^'^S,qg\^^^^ 



^S,gq = 



2X £/2 



+ 



JO) 1+ Jl) +£/,(!) 
75,99 g +«5,99 + ^"S,9 



+ 

+a 



1^(0) /^(O) , ^(0) \ _ Q 
275,99 1,75,99 + 75,99; Po 

^ (1) ofl (1) , (0) (1) , (0) (1) . ( da\ 
2rs,gq - 2/3o a^gg + 7^ ay^^ + 7^ a^_^^ " " 

9^ 1 



JO) 1 1 

^S,gq f 2 



J) ^ 
'S,gq 

da 



(2.26) 
(2.27) 



(2.28) 



(2.29) 



(2.30) 



2\^/2 



+ 



9 



,2 \2 



2\ <^ 



167r2 / /i2 



JO) i , Jl) , ..(1) 
75,99 £ +'*5,99 + ^''S,9i 

1 ('JO) f + 1 JO) JO) _/3,J0) 11 

2 V W +2^^'»«^^'«s ^075,99 J ^2 



^/-^J^) 9/? «W ^J°) «W ^J°) «W /j. 

+''5,99j<i=r 



> - 



(2.31) 
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Notice that in the above we have suppressed the Melhn index m. The expressions have been 
written in such a way that the anomalous dimensions take their values in the MS scheme. 
Furthermore we have in lowest order the identity "fs'qq = 7n^s qq- "^^^ above form of the 
unrenormalized OME's A^^ij follows from the property that the renormalized OME's A^^ij 
satisfy the Callan Symanzik equations ( |2.13| ), ( |2.14 ). These equations can be solved order 



by order in perturbation theory which provides us with the expressions presented at the end 
of this section. We only have to show that the latter follow from the renormalization of the 
OME's in p3^)-(|23ll). 



The renormalization of the OME's proceeds as follows. First replace the bare coupling 
constant g by the renormalized one g{fi) = g. Up to order g^ it is sufficient to substitute in 
the above OME's 

5 = 5(l + Yf^/3o5.^), (2.32) 



where (3q is given by ( 2.16| ). Next one has to perform gauge constant renormalization. Notice 



that in the next section we will calculate the one-loop OME's in a general covariant gauge. 
The Feynman propagator in this gauge is given by 

where a is the gauge constant. The two-loop OME's are computed in the Feynman gauge so 
that we have put in eqs. ( |2.25| )-( ^.3lD a = 1. 

Since the quarks and gluons are massless one has to put the external momenta p of the Feyn- 
man graphs off-shell. This implies that the OME's are no longer S-matrix elements and they 
become gauge (a) dependent. Therefore we also have to perform gauge constant renormal- 
ization which proceeds as follows. Replace the bare gauge constant a by the renormalized 
one. 



a 



ZnOt- 



(2.34) 



In the covariant gauge one has the property = Za where Za is the gluon field renormal- 
ization constant. Hence Za is given by 

g2 1 

Za = l + T^Zrv-, (2.35) 



167r2 e' 



where Za is given in (2.19). After these two renormalizations the only ultraviolet divergences 
left in the OME's are removed by operator renormalization. Choosing the MS scheme the 
operator renormalization constants are given by (see ( |2.21| ), ( |2.22D ) 

^2 \ 



Zns, 



1 + 



g 



167r2 



1 



(0) 

:7ivs, 



qq 



+ 



S 



ZNS,qq + ZpS,qq, 



qq 



with 



Zps,. 



qq 



167r2 



St 



2 



^S,qq'~^S,gq 



11 (1) 

+ —Tps. 



2e 



qq 



(2.36) 
(2.37) 

(2.38) 
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+ 



Zs. 



9Q 



+ 



(ifc) 

\167r2J 

167r2 



1 (0) ■ 



(2.39) 



2s ^'19 



1 



(0) 



m 



(2.40) 



(0) MO) , (0) \ , « (0) I 1 , 1 (1 



2£ 



99 



\ 167r2 



+ 



( fi'^ V c2 [/I MO) \2 , 1 (0) (0) . ^ (0) 1 1 
1^ 16^ J 2 W^'^J + 2 ^^i^^^is + /^o 7kg9 \ 72 



(2.41) 



Notice that the anomalous dimensions 'y^\j {k =NS,S, I = 0,1) are gauge independent so 
that Zk^ij have to be gauge independent too. The renormaHzed operator matrix elements are 
derived from 



and 



ANS,qq{p'^,l^'^,9,a) = Z^^ggg{g^,e)ANS,qq{p'^,l^'^,9,Ol)\g^Zgg; a^Z^a (2-42) 

(2.43) 



As,ijip ,IJ' ,g,a) = (Zs )ii{g ,e)As,ij{p ,m ,g,a)\g^Zgg; a^z^a 
with the results 



ANS,qq = 1 + 



167r2 



^^NS,qq^^[-^] + 



1^^ 



^NS,qq 



+ 



167r2 



1/(0) \ 2 _ 1 ^ (0) 
o\^NSm) APO^NS,qq 



InM-^ 



, ; 1 (1) (1) ,1 (0) (1) 



1 / '^"jVS.gq 



da 



-,(1) 



(0) ,(1) 



d 



da 



(2.44) 



a=l 



167r2 



99' 

1^(0) ^(0) , 2 
8 75,9ff 75,99^11 



+ 



8 



+ 



+ 



+«PS,g9 ^S,qg'^S,gq 



9^ 



167r2 ) 



'S,qg 



^ ^ 75,95 I75,gg + ^S,gg ) 4 75,gff / 1 ^2 I 
+ \2^S,qg - P0as,qg + ^ls,qg(^S,gg + 2^S,qq''S,qg 

+ «5,9ff + ^^005,95 ^S,qq^S,qg ^S,qgf^S,gg 

d ^ 



+ a 



167r2 



1 (0) , / V 

2^S,gq^^[-j2)+^S,gq 



,(1) 



\ 875,ffg l,75,gg + ls,gg) 4/^0 75,gg J I I 
+ j 2 ^5.9<? ^0"5,g<? + 2 ^S,gq(^NS,qq + 3 ^S,ggO'S,g 



-jlS,ggO-S,gq 



^ (1) (0) ,(1) (0) ,(1) 

+ (^S,gq + 2/^00^,^^ - ls,gq''S,qq " 7s,ag0s,99 



Y. 



a=l 



+ 



5 



.2 \2 



167r2, 



1 / (0) \ 2 , (0) (0) 1 ^ (0) 1 , 2 / 
+ \2^S,gg - Poas^gg + 2^S,gg^S,gg + 2^S,gq^S,qg 



2^S,ggO'S,gg+ 2^^'91^'19 



(2.45) 



(2.46) 



(2.47) 
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- as,59 + ^mbs^gg - ls,gg^S,9g " ^5,a<?^S,5g 



(2.48) 



The above renormalized OME's satisfy the Callan Symanzik equations in ( ^.131) , ( ^.14 ) which 



proves that the ansatz for the unrenormahzed OME's in ( 2.25 )-( 2.31 ) is correct. This is also 
corroborated by the expressions for the operator renormahzation constants Z^^ij in ( p.36 )- 
( p.41 ) which after insertion in eqs. ( p. 231) , ( |2.24 ) provides us with the anomalous dimensions 
in p30D . 

The above renormahzation procedure was originally introduced by F.J. Dyson There 
exists an alternative possibility invented by Bogoliubov, Parasiuk, Hepp and Zimmermann 
(BPHZ) lH]. In the latter one renormalizes each Feynman graph independently using the 
counter-term method. These counter-terms appear in the effective lagrangian which is ex- 
pressed into the renormalized (coupling- and gauge-) constants, fields and operators. The 
BPHZ-method has been used in the literature ||T^- p!7[| t o derive the anomalous dimensions 
of the spin averaged operators O^^'"^"^ (i=NS,S) in (|2.l| ). The advantage of this method is 
that the gauge dependent terms given by a '^^k\j) ^ol in ( p.25[ )-( p.31 ) are automatically 
subtracted. We will come back to this method at the end of section 3. The reason for the 
algebraic exercise given above can be explained as follows. Since the lowest order coefficients 
^^^\^ and a^^]j- can be very easily determined from the one-loop OME's one immediately can 
predict the double pole terms in the unrenormahzed OME's ( 2.25| )-( [2.30| ). The coefficient of 
the single pole term can be also computed except for the second order anomalous dimensions 
.(1) 



Ik ij- equating the predicted form of the two-loop OME's in (2.25)-(2.30) to the explicitly 



computed result in the next section one immediately can infer the results for 7^ 



(1) 



Calculation of the order contribution to the spin splitting 
functions 



In this section we first give an outline of the procedure of the calculation of the OME's defined 
in ( |2.11| ). Then we present the analytical result for the OME's and extract from them the 
splitting functions (anomalous dimensions). 

The calculation of the OME's proceeds as follows. Using the operator vertices in Appendix 
A and applying the standard Feynman rules we have computed the connected Green's func- 
tions, which are given by the graphs in figs. 1-6, up to two-loop order. The latter also involves 
the calculation of the diagrams which contain the self energies of the quark and the gluon 
in the external legs. These diagrams are not explicitly drawn in the figures but are included 
in our calculation. The computation of the one-loop graphs has been done in the general 
covariant gauge because one has to renormalize the gauge constant a even if one chooses the 
Feynman gauge a = 1. The two-loop graphs have been calculated in the Feynman gauge 
which is sufficient to that order. The OME's are then obtained by multiplying the connected 
Green's function by the inverse of the external quark and gluon propagators. Since the ex- 
ternal momenta are put off shell only ultraviolet divergences appear in the OME's which are 
regularized by using the method of n-dimensional regularization. This implies that we have 
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to find a suitable prescription for the 75-matrix which appears in the quark operators Rk,q 
for k = NS (p. 51) and k = S (pB). Here we will adopt the reading point method as explained 



in |20|. One can also adopt the method of 't Hooft and Veltman |]2l[, which is equivalent 



to the one given by Breitenlohner and Maison ||2^ (see also p^). The disadvantage of the 
last method is that the non-singlet axial vector operator R^j^^ g^ts renormalized in 

spite of the fact that it is conserved. This has to be undone by introducing an additional 
renormalization constant ||2^. However for continuity this procedure has to be extended to 
higher spin non-singlet operators R^s q {rn > 1) otherwise the anomalous dimension 

of R^sg ^^^^ become unequal to the anomalous dimension of the spin averaged non-singlet 
operator 0]v5g (^^^ (13))- Notice that the same procedure has to be also carried out for 
some of the singlet operators Rg'g ( ^.61) . Using the reading point method [^] one can omit 
the additional renormalization constant. Anyhow we have checked that both methods lead 
to the same result. 

As has been already mentioned in section 2 the Feynman rules for the operator vertices in 
Appendix A have been derived multiplying the operators -R^V by the sources J^^-^^ = 
A^^ • • • A^^ with A^ = 0. To simplify further we can choose s = p, where s is the spin 
vector in (|2.4| ), without any loss of essential information. The operator matrix elements ^^j^ 
( p.25D -( ^l3l[ ) are then given by 

Af^,^{p^^,^g,e){Apr = \Tr{GT,,{p,A,f,\g,eh,^}, (3.49) 
with k = NS, S and i = q,g and 

A^^,gip^^,^g,e)iApr = ^ ^^.a.A V G^f;(p, A, 5, e)- (3.50) 

Here G^- stand for the unrenormalized Green's functions which are multiplied by the inverse 
of the external quark and gluon propagators. 

We will now give a short outline of the calculation of Let us first start with the 

non-singlet OME A^g The Green's function G^g which is determined by the one-loop 
graphs in fig. la,b and by the two-loop graphs in fig. 2, consists out of Feynman integrals 
where the numerators are given by a string of 7-matrices. One of the 7 - matrices represents 
the 75-matrix. The latter is then anticommuted with the other 7-matrices until it appears 
on the right side of the string next to the 75 in ( 3.491 ). Then we set 75 = 1 and simplify the 



trace by contracting over dummy Lorentz-indices. Finally we perform the trace in (3.49). In 
this way one obtains the identity 

i^5,,,(p^^^e) = (i^5,,,(p^/z^.)) . (3.51) 

'^^ V '^^ / spm— averaged 

without any additional renormalization constant. Notice that the calculation of AJ^g^^ (spin 
averaged OME) has been already done in the literature so that it will not be repeated here. 
Except for the non-singlet operator A'^^^^ the remaining spin OME's differ from their spin 
averaged analogues. Since we need the one-loop OME's as presented in fig. 1, for the renor- 
malization of the two-loop OME's given by figs. 2-6 we have to calculate the former ones 
up to the non-pole term ^ defined in eqs. ( |2.25| )-( ^.31| ). The one-loop terms bk^ij which 
are proportional to e = n — 4, do not play any role in the determination of the anomalous 
dimension and they will not be presented in this paper. 
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Starting with the one-loop contribution to A^gq (fig- Ic) we have to perform tensorial 
reduction of the tensor integrals appearing in G^g^. These tensor integrals arise because 
the integration momentum appears in the numerators of the integrand. Examples of 
such one- loop integrals are given by eqs. ( [B.l| )-( p!3|) . The result will be that G^g^ gets 

terms of the form e°'^^'^ Aa-px^a^Jis, ^'^^^^ ^alalxlfj, ^^^^^Px^ala^lp, where the Levi- 
Civita tensor £°^l^^'' originates from the two-gluon operator vertex in ( [A. 41 ). Performing the 
trace in ( 3.49 ) provides us with a second Levi-Civita tensor so that we have to contract over 
two and three dummy Lorentz-indices. The contraction has to be performed in 4 dimensions 
since the operator vertices have a unique meaning in 4 dimensions only. Next we discuss the 
calculation of the one-loop contribution to A^qg ( 3.5C| ) (fig. ld,e). To this OME we apply 
the reading point method and put the 75 on the right hand side of the trace from the 
start. In this way we reproduce the Adler-Bell-Jackiw anomaly which can be traced back to 
the triangular fermion loop in fig. Id. Notice that fig. le leads to a zero result because the 
external momentum p appears twice in the Levi-Civita tensor. The Green's function 
will then become proportional to e^'^^'^ i^^xp^j. The latter will be contracted with the Levi- 



Civita tensor in ( |3.5C1| ) where the contraction is performed in 4 dimensions. The one-loop 
graphs contributing to ^'sgg ^^^^ presented in figs. If, g. Because of the Levi-Civita tensor 
coming from the two-gluon operator vector in ( |A.4| ) G^'gg will, after tensorial reduction, 
become proportional to e^^'^" A\Pfj. Like in the case of ^5"^^ the contraction with the Levi- 
Civita tensor in ( 3.5C| ) has to take place in 4 dimensions. 

Before we proceed with the two-loop graphs we want to emphasize that first the tensorial 
reduction has to be made before one can perform the contraction between the two Levi- 
Civita tensors. Both operations do not commute and lead to different results for the OME's. 
This holds for the one as well as two-loop calculation. If one contracts the Levi-Civita 
tensors in n dimensions both operations commute. However then the Lorentz indices of the 
operator vertices in Appendix A have to be generalized to n dimensions which is a non-unique 
procedure. 

The calculation of the two-loop graphs in figs. 3-6 proceeds in an analogous way as in the 
one-loop case. However here there arise some extra complications. First of all we encounter 
the two- loop scalar Feynman integrals which have already been performed in |24| to calculate 
the spin averaged OME A^gg- To check these integrals and the tensorial reduction algorithm 
we have recalculated all spin averaged anomalous dimensions (splitting functions) and we 
found complete agreement with the results published in the literature [p^-[17|. 
The second complication shows up in the tensorial reductions of the two-loop tensor Feynman 
integrals where the numerator now reveals the presence of two integration momenta qi and 
q2 ■ A more detailed explanation of how the tensor integrals are reduced into scalar integrals 
is presented in Appendix B. The third complication arises because of the appearance of a 
trace of six 7-matrices out of which two are contracted with the integration momenta qi and 
q2- Such graphs (see e.g. fig. 3 and figs. 5.11) are calculated by the following procedure. 
First one performs tensorial reduction of the Feynman integrals as indicated in Appendix B. 
This will lead to an increase of the pairs of 7-matrices having the same Lorentz-index. Then 
one can eliminate these pairs using the standard rules for 7-algebra in n dimensions. This 
is possible without ever touching the 75 matrix because it is put at the right hand side of 
the string of 7-matrices. After this procedure one ends up with the expression Tr{(fif^^^^) 
which is uniquely defined (irrespective of the 75-scheme). The same holds for the other graphs 
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which do not contain fermion loops. The final result is that all Green's functions get 
the same form as observed for the one-loop case. Four dimensional contraction of the two 
Levi-Civita tensors yields the OME's j4[™j in ( 3. 491) , ( |3.5[l| ). Before finishing the technical 
part of this section we give a comment on the algebraic manipulation programs which are 
used to calculate The matrix elements (including the full tensorial reduction) were 

calculated using the package FeynCalc [25] which is written in Mathematica ||2^ . The two- 
loop scalar integrals were performed by using a program written in FORM |25] which was 
called in FeynCalc. 

If one performs the inverse Mellin transform of the OME's the results for the one-loop 
calculation can be summarized as follows (see eqs. (2.21)-(2.27)). First we have the lowest 
order splitting functions which are already known in the literature ||T2| , p^ ]. 



P 



(0) 



NS,qq 



3(0) 
S,qq 



Cf 



X 



A-Ax + 65(1 



3(0) 

S,qg 



(0) 



S,gq 



P. 



(0) 



S,99 



Tf [I6x - 8] , 
Cf [8-4x], 

Ca 



X 



22 

16x + —5{1 - x) 



Tf 



^S{l-x) 



(3.52) 

(3.53) 
(3.54) 

(3.55) 



where the colour factors in SU(N) are given by Cp = (iV^ - 1)/(27V), Ca = N and Tf = 1/2 
{N = 3 in QCD). The non-pole terms a^j^jj appearing in expressions ( 2.25 )-( 2.31 ) read as 
follows 



NS,qq 



^S,qg 



^S,qq 



Cf 



ln(l 



X 



+ 2(1 +x)ln(l 



„l + x2 

2 mx 



X 



-4 + 2x + (l-a) 2 



1 



1 — X 

Tf [(4-8x)(ln3; + ln(l-x))-4], 
Cf [(-4 + 2x)(lna; + ln(l-x)) + 2-4x], 
ln(l — x) 



+ 5(l-x)(7-4C(2)) 



(3.56) 

(3.57) 
(3.58) 



Ca 



a 



1 — X 
1 



4 + 8x In X 



1 



X 



+ (-4 + 8x) ln(l -x) + 

+ 

+ 2 + S{1 - X) (^ - 4C(2) + (1 - a) - 1(1 - a)2) 



(3.59) 



In the above expressions the distributions i-x^'' ) defined by 

'ln'=(l-x)~ 



dx 



X 



fix) 



dx 



X 



im-fii)). 



(3.60) 



Notice that in the general covariant gauge only a^^gg (k=NS,S) and a.^s\g depend on the gauge 
parameter a. 
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The two-loop contributions to the unrenormahzed OME's are given by the inverse MeUin 
transforms (see eqs. ( p.25D -( p.31D ) 



^PS,qq 
^S,gq 

8 



(167r2)2 


St \ 








^-p2 


(167r2)2 










(16^2)2 




9' 







(167r 



2^2 e 



CpTfBf^, 

T] Bff + Cf Tf Bff + Ca Tf Bf^ 
Ca Cf B% + Cf Tf Bf^ + Cl B% 
Ca Tf B% + Cf Tf Bff + C\ Bfj, 



I 4(1 + x)lnx + 10(1 -x) I 
8 f 

+ - j 4(1 +2;)Li2(l - x) + 3(1 + x)ln2x + 4(1 +3;)lnxln(l - x) 
+10 (1 - x) ln(l - x) + (7 - 5 x) Inx - 5 (1 - x) I , 



n<79 
^// 


64 1 




^Ff 


- M 






_^4 


r _ 



(3.61) 
(3.62) 
(3.63) 
(3.64) 



- 8 (1 - 2 x) In X ln(l - x) + 3 (1 - 2 x) ln2 x + 2 (4 x + 3) ln(l - x) 
-(8x + 5)lnx-12 + 13x|, 

^ |-2 (1 - 2x) ln(l - x) + 4 (1 + x) Inx + ^(25 - 14x)| 
+^ I 12Li2(l - x) + 2(2x + l)Li2(-x) - 2 (1 - 4 x) C(2) 

- 3(1 - 2x) ln2(l - x) + 2(2x + l)lnxln(l + x) + 8 (1 + x) lnxln(l - x) 

2 1 11 

+(6x + 5) In^ X + -(11 X + 14) In X + -(73 - 62 x) ln(l - x) - -(44 - x) ^ 



^(25 X - 14) + 2 (2 - x) ln(l - x) - 2 (x + 4) In x| 
+^ I - 12 X Li2(l - x) - 2 (x + 2)Li2(-x) - 2 (4 - x) C(2) - 2 (x + 2) In x ln(l + x) 

+ 3 (2 - x) ln2(l -x)-{3x + 10) ln2 x - 2 (5x + 2) In x ln(l - x) 

1 111 
--(50 - 73x) ln(l -x) + - (17x + 8) Inx + -(109 - 119x) \ , 
o 3 9 J 
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T^gg 



^Ff 



^ |- (x - 2) j + - |- (10 - llx) + -{x - 2) ln(l - x) + - (x - 2)lnxj , 
^ I 3 X - 4 (x - 2) ln(l - x) + 2 (x - 2) In X I 

+^ I 12 (x - 2) Li2(l - x) + 8 (2 - x) C(2) + 6 (2 - x) - x) + 3 (x - 2) In^ x 
+4 (x - 2)lnxln(l - x) - 2 (4 - 7x)ln(l - x) + 7xlnx + 9 - 14x | , 

5(1 -x)| 

ln(l - x) + 3 



1 — X 



+ 2x - 1 



8 I 8 /ln(l -x) 
e j~3 V 1 -a; 

8 

+ 6x 



3 

2 ) Inx 



3 

26 



1 — X 



20 

3^ + y + 



/271 



C(2)U(l-x) 



3(1 - X 

^ I 4(1 + x) Inx + 10(1 -x) I 

+ ^ I 8(1 +x)Li2(l - x) +8(1 +x) lnxln(l - x) + 6 (1 + x) In^ x 

+20 (1 - x) ln(l - x) + (6 - 14x) In X + 22 (x - 1) + 5{l - x) | , 
'ln(l — x) 



16 



8 3 + 



1 — X 
1 

1 — X 



+ 16(l-2x) ln(l-x) + 22 



Inx + 20x - 42 + 5{l - x] 



( 



121 



1 

1 — X 
-8C(2) 



+- (-32 (1 + x)Li2(l - x) + f-16x - — ^ 8^ Li2(-x) 

el V 1 + X / 



+ 



+ 24 



-16x + 4 

/ ln^(l - X 
\ 1 -X 



1 — X 



+ 



+ 1 + 2;, 

88 /In(l-x) 



C(2) + -16x 



+ -48 X + 



+ 



8 



1 — X 



3 V 1 -X 
24 ) lnxln(l - x) + 



1 + x 

+ (-48x + 24) ln2(l-x 
2 10 



8 lnxln(l + x) 



/208 320\ 



119 29 



x) + -14x + 



88 



+^ - — X + 5(1 - x) — C(2) + 10C(3) 



1 + X 
88 
3(1 - x) 

1663 \ 



32 



1 — X 
38 ) Inx -43 



In^ X 



1 — X 



27 



Here the function Li2(y) stands for the dilogarithm which can be found in |25]. After sub- 
stitution of the one-loop order coefficients /3o ( |2.15| ), ( |2.19D , the Melhn transforms of P^Jj 
{= 'y^lf^^ ( 3.52 )~( 3.55 ) and a^j^jj ( 3.56| )-( |3.59| ) into the algebraic expressions for A^^ij in 

( p.25D -( ^.3lD one can equate the latter with the results obtained for A^^ij as presented above 
in (3.61)-( pj.64| ). From this one infers the two-loop contribution to the anomalous dimensions 
which are the unknown coefficients in eqs. (|]2|)-(|]3|). After performing the inverse Mellin 
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transform we get the sphtting functions. The non-singlet spHtting function Pj^l is the 



same as obtained in the spin-averaged case (see e.g. [|16|, |17[). In order to obtain the singlet 



splitting function Pg^gg one has to add to P^g the quantity Ppg given below. 



P\ 



(1) 



PSm 



Cf Tf -16(1 + x) In^ X - 16(1 - 3x) In x + 16(1 - x) 



(3.65) 



Furthermore we have the singlet splitting functions 



P, 



(1) 



ACaTj [-8(1 + 2x)U2{-x) - 8C(2) - 8(1 + 2x)lnxln(l + x) 
+ 4(1 - 2x) -x)- 4(1 + 2x) In^ x 
-16(1 - x)ln(l - x) +4(1 -F8x)ln2; - 44x + 48] 

+ACf Tf [s (1 - 2x)C(2) - 4 (1 - 2x) - x) 

+ 8(1 - 2x) lnxln(l - x) - 2 (1 - 2x) In^ x 

+16 (1 - x) ln(l - x) - 2(1 - 16x) In x + 4 + 6x] , 



(3.66) 



CaCf 16 (2 + x)Li2(-x) + 16 X C(2) + 8 (2 - x) - x) 
+ 16 (2 + x) In X ln(l + x) + 8 (2 + x) In^ x 

+ 16 (x - 2) In X ln(l ~ ^) + ( y + 3^ ) ^^(1 ~ ^) 

328 280 ■ 

+8 (4-13x)lnxH \ x 

9 9 

(x - 2) W{1 - x) - 4 (x - 2) In^ X - 164 + 128 x 
-8(x + 2)ln(l - x) - 4(20 + 7x)lnx] 



+cl 



+CFTf 



32 32 

■ — (4 + x) + y(x-2)ln(l-x) 



(3.67) 
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+CATf 
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-fsil-x) 



388 148 ... s 64, 
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(3.68) 
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For practical purposes, and the discussion of the results obtained above in the next section, 
it is also useful to present the one- and two-loop anomalous dimensions which are related to 
the splitting functions via the Mellin transform 



The one-loop contribution to the anomalous dimensions become 
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3 ^ 



(3.69) 

(3.70) 
(3.71) 
(3.72) 
(3.73) 



The two-loop non-singlet anomalous dimension 7]v5qq is the same as found for the spin 



averaged operator (see e.g. |15, |16[]). To obtain the singlet anomalous dimension jg'l'I^ one 



has to add the quantity Tp^'^ which reads 
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To check our results for the two-loop splitting functions (anomalous dimensions) we have 
also used the BPHZ method |]l9| as mentioned at the end of section 2. Here we renormalized 
the OME's graph by graph and found finally the same results as listed in ( p.65D -( p.67| ). As 
already mentioned in the beginning the above splitting functions and anomalous dimensions 
have been calculated in the MS scheme. If one prefers another scheme the corresponding 
anomalous dimensions are related to the MS ones in the following way 
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where ^k,qq {k = NS,S) denotes the anomalous dimension in the MS scheme and Z^s, Zij 
are finite operator renormalization constants. Up to order they can be expressed as follows 



Zns 
Z 



1 + 



9 



9 ^ 



^ + m^^99 



Substitution of eqs. (|3^ ), into eqs. (|]7|), dSlTOD yields 
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Ts. 



99^19- 



(3.80) 
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Before finishing this section we want to make a comment on the spin splitting functions 



and the anomalous dimensions calculated above. Two of them, i.e., P, 
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(1) ^.,(1) 



(7p5,gg) and 
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s qg (^Sqg) ^ave been already calculated in the literature They were obtained via 
mass factorization of the partonic cross sectoin of the subprocesses 'y* + q^q + q + q and 



gg 



'y*+g—>-g + q + q including the virtual corrections to ^* + g q + q. The result for Ppg 

( ^.65 ) agrees with eq. (3.37) in [^]. However the expression for P^^gg in ( 3.66 ) differs from 
the one obtained in eq. (3.38) of P] by a finite renormalization, i.e., 
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This finite renormalization is due to a different 75-prescription used in Q. 

The above splitting functions, which are calculated in the MS scheme, have to be combined 
with the quark and gluon coefficient functions ( 2.1C| ) computed in the same scheme in order 
to perform a complete next-to-leading order analysis. The quark coefficient function can be 
found in lll| and it equals to 
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The gluon coefficient function (see e.g. |2^) gets tlie form 
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The MelUn transforms of Ej. (k = q, g) become 



1 + 



167r2 



Cf 



+ - + 



m m + 1 



2 \ O"^ 

4Si{m-l)j In^ 

+ 3) 5i(m- 1) + 4 5i,i(m- 1) 



m m + 1 



6 



-AS2(m - 1) + 9 

m 



9 rji 

16^ ^ 



+ 



8 



m + 1 m 



m 771+1 



5i(m- 1) 



m m + 1 



(3.92) 



(3.93) 



Notice that the first moment E^ = 1 - 3{g'^ /{16Tr'^))CF agrees with eq. 6 of [10]. Further- 
more we have Eg = (see [0]). Both properties are characteristic of our choice of the 
75-prescription and the fact that the anomalous dimensions are calculated in the MS scheme. 



4 Properties of the spin anomalous dimensions 

In this section we will discuss some of the properties of the splitting functions and anomalous 
dimensions which have been calculated in the last section. Let us start with the first moments 
of the spin anomalous dimensions in the MS scheme. 

(4.94) 
(4.95) 
(4.96) 

^pTf, (4.97) 
(4.98) 
(4.99) 

where Pq and Pi are the first and second order coefficients in the perturbation series of the 
/3-function ( p^ ). 

In the above we have assumed that there is one flavour only in the fermion loops of the OME 
graphs. If there are more light flavours the Tj in the above expressions have to be multiplied 
by the number of light flavours indicated by (see (|2l^ ), (I2.17D ). The vanishing of the first 
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moment of the non-singlet anomalous dimension follows from the conservation of the axial 
vector current R^Sgq^^)- "^^^ value of the singlet anomalous dimension 'Js'^g was already 



calculated in [29|. It is due to the anomaly of the singlet axial vector current „{x) which 



contributes via the triangular fermion loops to Js'l'} in second order perturbation theory. The 



vanishing of 75 was shown on general grounds in [30|, see also |31]. From the last reference 



we also infer (see eq. (22) in |31]) that T^^^g = — 2/3i , provided the anomalous dimension is 
calculated in the MS scheme. Finally we want to investigate an interesting relation which is 
conjectured for an = 1 supersymmetric Yang-Mills field theory. It can be derived from 
QCD by putting the colour factors Cp = Ca = -/V and Tf = N/2 [32|. The relation reads as 
follows. First define 

= ISm + lS,gq - lS,qg " lS,gg- (4.100) 

For an iV = 1 supersymmetric Yang-Mills field theory one has 



57 = 0, 



(4.101) 



provided 5^ is calculated in a renormalization scheme which preserves the supersymmetric 
Ward identities. In many cases one has shown that at least up to two loops n-dimensional 
reduction is a regularization method which respects the supersymmetric Ward identities. 
Therefore a renormalization scheme where the pole terms plus the additional constants 
(Euler constant) and ln47r are subtracted (MS scheme) will respect these Ward identities 
too. In lowest order, where there is no difference between n-dimensional regularization and 
n-dimensional reduction, the above relation holds for the spin as well as spin averaged anoma- 
lous dimensions. If one assumes that the two-loop anomalous dimensions calculated in the 
two regularization schemes (n-dimensional reduction and n-dimensional regularization) are 
related to each other via a finite renormalization one can derive the following relation pF 



r (l),m 



(2 A 
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^S,qg 
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^S,gq 



'RED 
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(l),m 
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(4.102) 



with 



5a« 



^Sm + ^S,gq 



a 



Sm ^S,gg 



(4.103) 



where the terms a^'lj ^ieg ^^'^ "^s !j red non-pole parts of the OME's in ( ^.251) - 

( g3lD which are calculated using n-dimensional regularization and n-dimensional reduction 
respectively. Equation (|4.102| ) can be easily derived from eqs. (|!8|)-(|]8|) by putting ^7(0) = 



and z. 



V 



as,ij 



,REG^fls,ij,RED- If One makes the additional assumption ^7^^^ = ^7 



(1) 

RED 



( |4.10l| ) relation ( [4.102| ) turns out to be valid for the two-loop spin averaged anomalous 



dimensions which is checked in |33]. 

In the case of the spin anomalous dimensions we obtain the following results 
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(4.104) 
(4.105) 
(4.106) 

(4.107) 
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If we assume that Sjj^^^y) = then (from ( 4.102D ) 



(1).™ 



+ 



28 1 44 1 

+ 



(m + 1)2 3 m 3 m+l' 



(4.108) 



which is in disagreement with the result of our calculation derived from eqs. ( |3.74D ( |3.77| ) 
which is equal to 



c- (l),m , , 



548 1 532 1 



3 m 3 m + 1 



104 



(m + l)2' 



(4.109) 



We tried to explain the difference between the prediction in (|4.108|) , which is based on su- 
persymmetry, and the result ( [4.10S| ) obtained by our calculations. Therefore we investigated 
the findings in [33| for the spin averaged anomalous dimensions and found a surprising result. 



Since all external legs of the OME's are put off shell one can split G^^g ( |3.49| ) and G^'-g''^ 



( p.50D into a so-called physical and unphysical part. In the case of G^jg the former part is 
proportional to Z)( whereas the latter part is multiplied by ^. This property holds for the 
spin as well as spin averaged operators. In the spin case G^fg^ has a physical part only 
which is proportional to £^^"l^ /S.aPi3- However for the spin averaged case one also encoun- 
ters unphysical parts in the OME's. Here the physical part is the coefficient of the tensor 
— {IS.^py +p^Ajy)/(Ap) + A^Aj^p2/(Ap)2. Limiting ourselves to the physical parts of 

the non-pole terms CLg\j we find the following results. The spin averaged OME's satisfy the 
relation (Feynman gauge) 
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whereas the spin OME's lead to 
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RED 
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(4.110) 



(4.111) 



Property ( 4.110 ) was not mentioned in |33]. However it might explain why = for 



the spin averaged case since the physical part of the unrenormalized one-loop OME's already 



satisfy the supersymmetric relation. Hence we have the suspicion that because of ( 4.111 ) 



"^Tred 7^ ^ ^P™ anomalous dimensions which explains the discrepancy between ( [4.108| ) 



and ( [4.109 ). Notice that ( [4.11l| ) is obtained in the MS scheme. Therefore we have made an 
oversubtraction so that ^Oj^gj-, = 0. If we now assume that after this subtraction 5'^'^y) — 



oversubtraction so that 5a 

and recalculate 5a^^Q (Feynman gauge) by keeping the physical part only we obtain 
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(4.112) 



and from ( [4.102| ) 



•^Treg 
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+ 



24 100 1 116 1 

(m + 1)2 3 m 3 m + 1 



(4.113) 



By comparing ( 4.10S| ) with ( 4.113| ) we observe again a discrepancy between our calculation 
and the prediction obtained from the supersymmetric relation except for m = 1 where we 
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have agreement. The reason for the violation of the supersymmetric relation (4.101) in the 
case of the spin anomalous dimensions is not known to us. It cannot be attributed to an 
error in the scalar Feynman integrals because they were also used to recalculate the spin 
averaged anomalous dimensions which we found to be in agreement with the results quoted 
in the literature. It might be due to our 75-prescription. However different 75-prescriptions 
are related via finite renormalizations and the discrepancy between ( [1.1 09| ) and ( 4.108| ) or 
( [1.1 13| ) cannot be explained by such an effect. Finally we want to emphasize that to our 
knowledge the formal proof of the supersymmetric relation ^7 = is still lacking in the 
literature. 



Appendix A: The operator vertices 



In this appendix we present the twist-2 operator vertices. All momenta are flowing into the 
operator vertex. 



A.l Quark- (gluon) operator vertices 

The quark-antiquark vertex is equal to 

o{p) = -4j,{Apr-\ 



(A.l) 



where p denotes the momentum of the incoming quark line. 
The quark-quark-gluon vertex is given by 



m—2 



Oil{p,q) 



A^4(75 5: (Ap)— 2(_Agy 



(A.2) 



i=0 



where p and q are the momenta of the incoming quark and antiquark respectively. 
The quark-quark-gluon-gluon vertex equals 

O^J{p,q,r,s) = g^A^^A-^/^j, 

(m-3 j 
(-1)" E E {-lyi^pYiMr-'-M^p) + (my-' 

/m-3 j \" 

-TbTa E E i-iyi^pr-'-HMynAq) + (a.))^- . (a.s) 

\j=Q i=0 )_ 

A.2 Gluon-operator vertices 

The 2-gluon vertex is given by 

OZip) = te^^'^^il - {-ir){^pr-^5ab. (A.4) 
The 3-gluon vertex is equal to 
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with 

—[e^p^'^A'^ — eP-P^"^ A^){Ar)'^~'^ 

m— 3 

+ ^ {-Apy{Aq)'^-'-^AP{e''P^iA^ + e^"'^1{Ap)) 

i=0 
m-3 

- J2 {-Ary{Ap)"'-'-^A''{eP-P^'"AP + e^P^PiAr)) 

1=0 
m-3 

- J2 {-Ary{Aq)"'-'~^Ai'{e''^'i''AP - e''P^'i{Ar)). (A.6) 

The 4-gluon vertex equals 

0'aZ{p,q,r,s) = igHl-i-ir)[fabefcdeOP''P''ip,q,r,s) 

+facehd.OPP'"' {p, r,q,s)- fadefbceOP'^P'' (r, q,p,s)], ( A.7) 

where 

0^''P''{p, q, r, s) = {e^'^P^A^' - e^^P" A^'XiAr) + (As))"^-^ 

m— 3 

_^p(^^uaAs^l, _ ^,xaAs^u^ ^ ^ (As))^(As)™-^-^ 

i=0 
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^^a^^pz/r-A^M _ gP^^'^A^i.) ^ ((-^^^ ^ (A9))™-^-3(-(Ar))^ 

i=0 
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+A'^ ^ ((Ar) + {As)r-'-^i-iAp)y{e'"'^PAP - e'^p^pA'') 

i=0 
m—3 

+A'^ J2 ((Ar) + (As))"-*-3(-(A9))*(e^^^«A'' - e'^P^iA") 

i=0 

m-4 i 

+AW ^ ^(Ap)™-^-^((Ap) + (Ag)y-*(-(As))*(e^"P^A'^ + eP''^'{Ap)) 

m-4 j 

-A^AP J2 J2(^ir~^~^((^P) + (Ag)y"*(-(As))*(e^"«*A" + e'^"^^(Ag)) 

j=0 i=0 

m— 4 J 

-A'^A" ^ ^(Ap)"^-J-^(-(Ar))*((A;j) + (Ag))^-X£^^^'"A'' + e^''^P(Ar)) 

j=0 j=0 

m— 4 j 

+A'^A<^ ^ ^(Ag)"^-^-4(-(Ar))*((Ap) + {Aq)y-\e'^'"i'^AP + e^'P^iiAr)). (A.8) 
Appendix B: The tensorial reduction 

In this appendix we present a more detailed explanation of the tensorial reduction of the 
tensor Feynman integrals into scalar integrals. 
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According to the reading point method we can put the 7^-matrix at the right end of 
the traces. Then one can perform all straightforward simplifications of the 7-matrix algebra 
inside the traces. Furthermore we leave the 7^-matrix untouched except that at the end 
we take Tr ^7°^ 7^ 7"^ = —4 i e"^°"^ In the case of the one- loop integrals the tensorial 
reduction can be very easily achieved via the standard Feynman parameter techniques. Since 
we do not need tensor integrals beyond rank two it is sufficient to list the following integrals 
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(B.3) 



where 5„ is the spherical factor 5„ = vra /(27r)". 

The tensorial reduction of the two-loop tensor Feynman integrals is much more complicated 
and has been performed by using the program FeynCalc |31]. The numerators of the two-loop 
Feynman integrals have the following structure. 

A = fr''''ha,-.Muq2). (B.4) 
where qi and q2 denote the integration momenta. Explicit forms of fp'"^'' and /j,o-i ■••0-^(91' Q2) 
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are (e°^'^VA/3 = eP^''^ etc.) 

/i (91,92) = e''^^^^^ /i(gi,92) = e'^^'^^'^P, (B.5) 

/2"(ft,92) = e"^^^", ka{qi,q2) = e"^^^", (B.6) 

/3"^(9l,g2) = e'?^^^"^ /3,a/3(ft,'Z2) = TrWi^24V7a7/375). (B.7) 



The tensor integrals can be represented as 



{«"; of 9?; ifv^j€:'i'v^j€4} , (B. 



where i = 1, 2 and j = 1,2. Further we have the definition 

J., . . _ (Agi)"(Ag2)''(A(p - gi))^(A(p - g2))'^(A(gi - 92))^ 

Notice that the integers a — g can take positive as well as negative integer values. By virtue 
of Lorentz covariance the integral I{A,p) can now be written as 



(27r)" 7 (2^ 



X 



I^I ^Efrip',n)Mquq2,A,p), (B.IO) 



with 



(B.ll) 

'A^}, (B.12) 
"A'^A'^}, (B.13) 
i°A^A'"A^} , (B.14) 

where the Kr are of the same type as the K in ( p. 8] ) (but with different indices a-j) and the 
fr{jp^,n) are simple polynomial-like functions determined by the tensorial reduction. 
In this way all Lorentz indices are transformed away from the integration momenta to the 
external momentum p and the lightlike vector A. The advantage of the tensorial reduction 
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method is revealed when one evaluates e.g. the expression fs^ap (B.7). This gets simplified 
to Tr(Z)(^ 7^7/375) = — 4ie^^"^. Hence one can avoid any 75-prescription dependence arising 
from the non-unique way of calculating a trace of six 7-matrices plus the 75-matrix in n 
dimensions. The explicit reduction formalae, which are too lengthy to be presented here, are 
obtained by using projection methods. They are incorporated in the program FeynCalc 3.0 



1 25]. The scalar integrals which appear on the right hand side of ( B.IO ) are calculated in [10| 



using the algebraic manipulation program FORM |27]. The two-loop integrals including the 
tensorial reduction have been checked by recalculating the spin averaged splitting functions 
which have been computed in the past (see [|l4|-Il^]) and we found full agreement. 
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Figure captions. 



Fig. 1 One-loop graphs contributing to the spin OME's; (a), (b): ^7^5 gg, A^s\q'-, (c): ^5!,^; 

(d), (e): A^s\g'i (^)' (s)- Graphs with external self-energies and with triangular 

fermion-loops where the arrows are reversed have been included in the calculation but 
are not shown in the figure. Graphs which arc not symmetric with respect to the vertical 
line through the operator vertex have to be counted twice. 

(2) 

Fig. 2 Two-loop graphs contributing to the spin non-sing let OME Graphs with 

external sclf-encrgics have been included in the calculation but are not drawn in the 
figure. Graphs which are not symmetric with respect to the vertical line through the 
operator vertex have to be counted twice. 

(2) 

Fig. 3 Two-loop graphs contributing to the spin pure-singlet OME Apg^^. Graphs with 
triangular fermion loops where the arrows are reversed have been included in the cal- 
culation but are not shown in the figure. Graphs which are not symmetric with respect 
to the vertical line through the operator vertex have to be counted twice. 

(2) 

Fig. 4 Two-loop graphs contributing to the spin singlet OME Ag'^^. Graphs with triangular 
fermion loops where the arrows are reversed and diagrams containing external self 
energies have been included but are not shown in the figure. Graphs which are not 
symmetric with respect to the vertical line through the operator vertex have to be 
counted twice. 

(2) 

Fig. 5 Two- loop graphs contributing to the spin singlet OME ^5^^- Graphs with external 
self-energies have been included in the calculation but are not drawn in the figure. 
Graphs which are not symmetric with respect to the vertical line through the operator 
vertex have to be counted twice. 

(2) 

Fig. 6 Two-loop graphs contributing to the spin singlet OME ^5^^- Graphs with external 
self-energies and diagrams with ghost and triangular fermion loops where the arrows are 
reversed have been included in the calculation but are not drawn in the figure. Graphs 
which are not symmetric with respect to the vertical line through the operator vertex 
have to be counted twice. 
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Figure 3: 



Figure 4: 
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Figure 6: 
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